Quantum light generation on a silicon chip using waveguides and resonators 
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Integrated optical devices may replace bulk crystal or fiber based assemblies with a more compact and con- 
trollable photon pair and heralded single photon source and generate quantum light at telecommunications 
wavelengths. Here, we propose that a periodic waveguide consisting of a sequence of optical resonators may 
outperform conventional waveguides or single resonators and generate more than 1 Giga-pairs per second from 
a sub-millimeter-long room- temperature silicon device, pumped with only about 10 milliwatts of optical power. 
Furthermore, the spectral properties of such devices provide novel opportunities of wavelength-division multi- 
plexed chip- scale quantum light sources. 



Trends in quantum optics are evolving towards chip-scale 
photonics [1], with one of the eventual goals being the full- 
fledged combination of sources, circuits, and detectors on a 
single chip. Regarding chip- scale sources, researchers have 
predicted and shown that an optically-pumped spontaneous 
four-wave mixing (SFWM) process in silicon can be used 
to generate entangled photon pairs in waveguides and res- 
onators J2-[l] • This third-order nonlinear process is similar to 
the second-order spontaneous nonlinear optical processes in- 
duced in bulk optical crystals (except scaling with the square 
of the pump power instead of linearly), and before being in- 
vestigated in lithographically-fabricated waveguides, has been 
demonstrated in optical fiber 10, 0]. As a further step, we have 
explicitly shown heralded single photons at 1.55 jum wave- 
length from a silicon chip at room temperature (see Fig. [T|) [8]. 
Given the maturity of integrated optics technology, it is realis- 
tic to envision on-chip high-brightness single-photon sources 
at wavelengths compatible with the worldwide fiber optic in- 
ternet infrastructure. 

However, an important open question is: What is the opti- 
mal device for generating quantum light using an integrated 
photonic structure? To be specific, we focus on devices 
made using silicon. The photon pair generation rate de- 
pends on the intrinsic four- wave mixing nonlinear coefficient 
(y = 271^2 /XA e ff), in terms of the Kerr nonlinear index ri2 and 
the effective area of the waveguide mode (A e ff), the waveguide 
length (L), the pump power (P), and the loss coefficient of 
lithographically-fabricated waveguides (a). Silicon nanopho- 
tonic waveguides are already quite promising, compared to 
optical fiber or bulk crystals, since a single mode " wire" 
waveguide with cross-sectional dimensions of about 0.5 x 
0.25 jum 2 has a nonlinearity coefficient ysi= 100-200 W _1 m _1 
(five orders of magnitude greater than optical fiber) around a 
wavelength of 1.5 jum J9|]. But chip-scale devices present spe- 
cial challenges as L is limited to only a few centimeters within 
a typical die site, and on-chip footprint is a highly valuable re- 
source in CMOS fabrication. Moreover, for a waveguide that 
is fabricated with loss coefficent a, the effective interaction 
length of nonlinear interactions L e fj = [1 — exp(— aL))/a can 



be significantly smaller than L when aL > 1 . Also, pump pow- 
ers P in silicon are limited to a few milliwatts to minimize the 
probability of multi-photon generation and avoid two-photon 
absorption and free-carrier generation losses. 

The indistinguishability of output single photons is also 
an important consideration [10]. In silicon waveguides, the 
phase-matching bandwidth of the SFWM process is generally 
quite broad, on the order of tens of nanometers. As such, the 
generated photon pair emerges as an entangled state, and de- 
tection of the heralding photon projects the signal photon into 
a mixed state. Purity may be enhanced by spectrally filter- 
ing the output, the disadvantage being a reduction in photon 
count rate since unused pairs are discarded. Recent work ifTTIl 
has also shown that through the careful control of waveguide 
dispersion, photon pairs may be generated in factorable states 
which are spectrally de-correlated. Alternatively, one may 
limit the modes available for SFWM process by placing the 
nonlinear material in a cavity, thereby providing both spectral 
filtering of output states as well as local intensity enhancement 
of the pumpfijQ. 

Based on these considerations, we study whether a micro- 
resonator is preferable to a conventional waveguide as a her- 
alded single photon source with specific reference to silicon 
devices. We then show that a particular type of hybrid device 
[Fig. [TJ, which consists of a linear array of nearest-neighbor 
coupled microresonators, can possibly generate in excess of 1 
Giga-pairs per second for 10-20 mW of optical pump power, 
from a waveguide that is only 0.1 mm long, thus outperform- 
ing existing photon pair sources by 1-2 orders of magnitude 
in generation rates and by 2-3 orders of magnitude in device 
size. 

In single ring resonators, the theory of both parametric 
downconversion (second order nonlinearity) as well as SFWM 
(third order nonlinearity) has been studied iTTsI — TTvh . Here we 
extend the previously described methods to develop the output 
state of the photon pair from a series of directly coupled rings, 
so that waveguides, rings and coupled-ring waveguides can be 
compared. We begin with the phenomenological Hamiltonian, 
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where aj m are the field operators of the resonator modes 
/ = s, i at the resonator site m, £l m are the resonance frequen- 
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Figure 1: (a) A coupled resonator waveguide consisting of N 
directly-coupled microring resonators. The waveguide eigenmode 
is a Bloch excitation, i.e., a collective oscillation of all N resonators, 
with a fixed relationship between adjacent resonators II 1411. The di- 
rection of light circulation in each resonator is as indicated for the 
specified input. In the notation used in this paper, the field operator 
of successive resonators are a\, ai, 03, the resonance radial fre- 
quencies are £22> £^3> the inter-ring coupling coefficients are 
labeled K2, K3, K4, and the input/output external coupling coeffi- 
cients are labeled by their rates J- and J- (the latter is not shown, 
at the output side of the chip). The resonator loss is indicated by 
the damping rate ^. (b) Light from a laser diode (LD) is coupled 
into the waveguide device under test (DUT), and output photon pairs 
are spectrally separated (wavelength division multiplexer, WDM), 
filtered (band pass filter, BPF), and detected e.g., by single photon 
avalanche detectors (SPADs) shown here in the configuration needed 
to perform a (0) measurement using a time-correlated single pho- 
ton counting system (TCSPC) and a data acquisition card (DAQ). 



cies, K/ ?m+ i are the inter-resonator coupling coefficients and 
% m is the coefficient proportional to the Kerr nonlinearity. In 
general % m may be time-dependent, % m (t) = ^p- [A p ^ m (t)] 2 , 



where A p ^ m (t) = a PiYn (t)e iQ, p* is the classical pump with a 



slowly- varying amplitude at carrier- frequency £l p , Jo is the 
usual waveguide nonlinear parameter, v g is the waveguide 
group velocity, T c = 1 / FSR is the round-trip time (inverse of 
the free-spectral range). We adopt the approach of Collett and 
Gardiner lll8ll (i.e., time-domain coupled mode theory) to ob- 
tain the equations of motion in the Heisenberg picture. In the 
single resonator case, these may be written explicitly as, 
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where a s ((O s ) are the frequency components of the time- 
dependent field operator a s (t) and ^ = ^ + ^ is the damping 
coefficient which includes effects of loss and external cou- 
pling. These equations contain the same information as the 
joint- spectral amplitude, modified by the cavity enhancement 
effects. In the quasi-cw limit, one may forgo the integral and 
solve the coupled equations as was done in II 1911 . 
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We have used the boundary condition \a out \ 2 = jn 2 \a\ 2 , where 
/j 2 = y is the input mode coupling coefficient. In the 
case of vacuum input and low gain the power spectral den- 
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and the total signal flux is F 



2^ f o(co)dco, where the idler frequency is implicitly related 
by the energy conservation 2co p = CO, + co ; . Alternatively, by 
taking %(co, + CO;) as the pump distribution in the pulsed pump 
regime, o(co,,COz) is interpreted as the joint spectral intensity. 

Extending to the case of multiple coupled cavities |20|], we 
have the following matrix equation, 
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and we have assumed a single sided input/output. 

Similar to single ring case, we have for the coupled- 
resonator waveguide, 

aoutA^s) = -JUlH2[TN,l(tos,®i)ain,s(to S ) 

+ TV^+i (©„ co/)a^.((o,-)] (5a) 
+ ?2A^+i (©* , i(<*>i)] (5b) 



and the joint spectral intensity o(cOs,Cfy) 
We note here that the coupled mode theory re- 
sult is equivalent to the first-order perturbation the- 
ory with a cavity modified joint spectral amplitude, 
l¥> = \0) 5 \0)i + gffd® s d<Bim s (<Q s )FEi(<Qi)FE p ((i) 5 ,(i)i) x 
/((0 > s,co i -)^(co > s)^(Q) i -)|0) >s |0) i - where the subscripts s and i 
refer to the signal and idler frequencies, g is proportional to 
the photon-pair production rate, and the function /((Q 5 ,(Q;) 
which describes the phase-matching and pump spectral 
envelope, is the joint spectral amplitude lfl3ll . FE are the field 
enhancement factors and are equivalent to the slowing factors 
used in |2lll . 

We verify the agreement between the time-domain coupled 
mode equations and the slowing factor enhanced pair gener- 
ation equations by comparing the calculated pair flux. In the 
discussion below, we will assume a simplified picture with flat 
spectral filtering about the desired signal and idler modes, as 
was done in previous experiments Il22ll . The number of pho- 
ton pairs generated per second is given in the low pump power 
regime by 



F = Av (Yeff/>Leff) exp(-aL) 



(6) 



where = S s Sj -(J^j-^ 7q> S{p,s,i} are me slowing factors 
at the pump (/?), signal (s) and idler (i) wavelengths, and 
L&s = [1 — exp(— aL)]/a represents an effective propagation 
length, defined as the geometric length L = N%R renormal- 
ized by the loss coefficient, a. An experimentally-validated 
transfer-matrix method can be used to calculate the a co- 
efficient which scales linearly with the slowing factor l23ll . 
We assume that the linear loss coefficient a does not vary 
significantly with wavelength over the bandwidth of inter- 
est. To account for nonlinear absorption losses in silicon 



J!] we substitute a a + 2^(3L and PL Q& PL where 

T = [log(l + ;Jff PL eff)]/l|f L and P is the effective TPA co- 
efficient of the coupled-resonator waveguide which scales in 
the same way as y e ff with S, i.e. (3 oc S 2 $o. For an apodized 
structure, which we define as the case where the boundary 
coupling coefficients are matched to the input/output waveg- 
uides [24], we have at resonance S = 1/|k|, where |k| is the 
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Figure 2: (a) Calculated photon pair flux F using pair generation 
equations, Eq. J6]). The white trend- line follows the optimum num- 
ber of resonators for a given slowing factor, (b) Corresponding values 
°f JeftPL for each S and N, showing the low multiphoton generation 
probability along the white line, (c) Calculated photon pair flux F us- 
ing coupled mode equations, Eq. flat . The top region of the contour 
plot represents a single resonator, while the far left approaches that of 
a conventional silicon nanowire waveguide. For S = 50, the optimum 
number of resonators is N opt = 25 for which F = 4 MHz/mW 2 . 
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inter-resonator coupling coefficient in the transfer-matrix for- 
malism. The bandwidth of the photon generation process, Av, 
is assumed to be the linewidth of a Bloch eigenmode of the 
coupled-resonator waveguide, which scales inversely with the 
number of resonators in the chain, N. 
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Calculations were performed using the following parame- 
ters, R = 5 jum, waveguide loss = 1 dB/cm, yo = 200 W _1 m _1 , 
Po = 0.75 cm/GW, P = 1 mW to obtain F over a range of val- 
ues of S and N, showing good agreement between the pair 
generation equations and coupled mode equations [Fig. [2] (a) 
and (c) respectively] . Resonator chains that are in excess of 
the optimum length, or with too high a value of S incur penal- 
ties because of the exponential loss factor in Eq. ©, and 
the collapse of the bandwidth Av. Too small values of S do 
not fully utilize the slow-light enhancement of the nonlinear 
FWM coefficient, which scales as a higher power of S than 
the corresponding decrease of bandwidth, unlike in a (linear) 
slow-light delay line. The optimum parameters are large S 
and small N, i.e. towards the single resonator configuration, 
for which the maximum pair flux rate exceeds 10 MHz at 1 
mW pump power (and scaling quadratically with the pump 
power). 

Scaling difference between single rings and coupled ring 
waveguides: One of the questions regarding the optimum de- 
vice geometry for generating photon pairs is the appropriate 
size of resonators. Recently, the efficiency of classical and 
spontaneous four- wave mixing in single microring resonators 
has been compared J25l|2^], with the conclusion being that in 
both cases, the conversion efficiency scales with the ring ra- 
dius as R~ 2 , i.e., smaller rings are better than larger rings in 
generating photon pairs. This results from the analytically de- 
rived expression for the spontaneously-generated idler power 
Pi,sp (from an injected pump power P p at optical carrier fre- 
quency co p ) 
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and a key assumption, that the ring quality factor Q is inde- 
pendent of the ring radius R. Starting with the equation for 
the (loaded) quality factor of a ring resonator side-coupled to 
a waveguide Il27ll . 
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where a rt — exp(— ocL/2) and L = 2nR, we examine two lim- 
iting cases as examples. In the first case, we examine a 
weakly coupled resonator (x = — |k| 2 = 1) with low loss 
(a n ~ 1 — ocL/2) in which case the quality factor can be ex- 
pressed as, 
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which is the intrinsic Q limit. In this case, Q is indeed in- 
dependent of R, and P^sp scales as R~ 2 . In the second case, 



however, we assume that the loaded Q is dominated by the 
coupling coefficient (|k| ^ 0) and then 
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In this case, the ratio Q/R in Eq. ([8]) is length-invariant, and 
P^sp increases linearly with R. As previously shown j28ll , 
coupled-resonator waveguides are more disorder tolerant in 
the large-coupling regime, and therefore, Eq. d9t) is more ap- 
propriate in describing performance rather than Eq. ^p). In 
fact, the agreement between Fig. [2a), calculated using the 
conventional waveguide theory with nonlinearities scaled by 
the slowing factor, and Fig. [2c), calculated using the first- 
principles time-domain coupled-mode theory model, shows 
that coupled-resonator waveguides are, in fact, more similar 
to waveguides than single resonators in many ways, with the 
attendant benefits of a slowing factor in enhancing the nonlin- 
earity per unit length. Here, it is useful to recall, as shown 
in the classical domain, that coupled-resonator waveguides 
break the traditional trade-offs between parametric conversion 
efficiency and bandwidth, and are more robust against chro- 
matic dispersion and propagation loss, compared to conven- 
tional waveguides [29] . Similarly, in the quantum domain, 
coupled-resonator waveguides may outperform conventional 
waveguides as pair and heralded single photon sources. 

Multi-photon generation probability: For a heralded single 
photon source we require low multi-photon probability. The 
inset of Fig. [2ta) shows the value of the quantity y e f[PL for 
each value of S and N. For a y e ff PL <C 1 , the level of stimulated 
scattering events is kept relatively low yfl which is true for the 
regions of highest pair flux (large S and small N). 

Joint Spectral Intensity (JSI): To evaluate the spectral char- 
acteristics of the signal-idler photon pair, we calculate the 
Joint Spectral Intensity, and also the Schimdt number K = 
1/ HX 2 , which is the sum of the squares of the Schmidt eigen- 
values (for a pure state K = 1) 113011 . In Fig. Oa), (b), we plot 
the joint spectral intensities of an unapodized and apodized 
coupled-resonator waveguide of similar inter-resonator cou- 
pling coefficients. The shape of the spectrum reflects the num- 
ber of resonators chosen N = 5, with the peaks corresponding 
to the locations of maximum transmission, which are also the 
Bloch eigenmodes. The pump pulse width is taken as 10 ps 
in both cases and we obtain K = 4.47 for the unapodized de- 
vice and K = 3.31 for the apodized device. However, we note 
that choosing shorter pulses does not significantly change the 
Schmidt number in contrast with the single ring case 11611 . In 
order to herald pure state single photons, filtering will be nec- 
essary. Choosing a filter bandwidth equal to the Bloch eigen- 
mode width given by Eq. ©, we are able to obtain approxi- 
mately a single Schimdt mode output. 

On the other hand, if we have control over each individ- 
ual inter-resonator coupling coefficients we are able to syn- 
thesize a large variety of different joint spectral amplitudes 
with different Schimdt numbers. In Fig. [2c), (d) we plot two 
interesting contours taken from a sample of different inter- 
resonator coupling configurations, each coefficient being a 
pseudo-random number ranging from to 1. Clearly, with 
the added control over individual couplers we can obtain a 
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Figure 3: Joint Spectral Intensity (JSI) plots for various coupling co- 
efficient configurations, assuming that the coupling coefficients be- 
tween adjacent resonators, shown in Fig. 1, can be individually al- 
tered, (a) Unapodized (b) Apodized (c),(d) Chosen from a sample 
of Monte Carlo simulations with random coupling coefficients, (e) 
JSI for coupling coefficients chosen so as to realize a Butterworth fil- 
ter response and (f) Bessel filter response in the linear transmission 
regime. 



large variety of corresponding K values. Of special inter- 
est are the configurations giving maximally flat transmission 
(Butterworth) and maximally flat group delay (Bessel) l24ll 
since these quantities define the overall shape of the output 
joint spectrum [see Fig. Efe), (f)]. Without additional filter- 
ing, we are able to obtain close to a pure heralded state for 
both the Butterworth filter configuration (K = 1.18) and the 
Bessel filter configuration (K = 1.09). Of course, filtering will 
still be required before the detectors, to separate the signal and 
idler photons and reject any unused pumps from reaching the 

spads (i. 

While we have focused on the details of a single resonance 
in the prior discussion, as was predicted for for the case of a 
single resonator lfl7ll . the full two-photon state generated by 
the coupled resonator device is expected to form a "comb" 
structure with peaks centered around the resonance frequen- 
cies. In Fig. Ua) we plot the transmission spectra around 
five particular resonances of a 5 -ring unapodized coupled res- 
onator waveguide, taking into account both the dispersion of 
the intrinsic constituent waveguides as well as the dispersion 
of the directional couplers Olfl . The spectrum of the two pho- 




Figure 4: (a) Spectra of the transmission bands of a coupled resonator 
waveguide consisting of five microrings. (b) Spectrum of the two 
photon state when a cw pump is placed at the resonance £l p . (vertical 
axes are in logarithmic scale for both (a) and (b)) (c),(d) JSI of the 
transmission bands adjacent to the pump as well as two bands away. 



ton state for a cw pump placed at the resonance £l p is given in 
Fig. HJb), showing a fine structure characteristic of the number 
of resonators. While the general structure remains consistent, 
the peaks near the edges are reduced more quickly than those 
near the middle. This can be attributed to the large directional 
coupler dispersion which give rise to non-uniform transmis- 
sion bandwidths. Careful inspection of Fig. HJa) shows that 
the bandwidths increase gradually with frequency. The fur- 
ther apart the bands are, the more misaligned the transmission 
peaks become which in turn reduces the effective nonlinearity 
(see Eq. [6]), since transmission peaks correspond also to peaks 
in slowing factor. The band edge peaks are most adversely af- 
fected since they are also the narrowest. In Fig. 0tc,d), we plot 
the JSI with signal and idler in the adjacent resonances as well 
as being two resonances apart from the pump. As compared 
to Fig. Efa), we can see that the band edge peaks have become 
more distorted. Clearly, the uniformity of the two photon state 
generated over the "comb" for the coupled resonator configu- 
ration is limited by the dispersion of the directional couplers, 
the suppression of which is a problem of interest not only for 
chip-scale quantum optics but in "classical" photonics as well. 

In summary, we have calculated the expected photon pair 
flux rates from a silicon coupled-microring waveguide de- 
vice based on spontaneous four- wave mixing, a nonlinear pro- 
cess which scales quadratically with the optical pump power. 
This hybrid structure may significantly outperform conven- 
tional waveguides of much longer length at realistic waveg- 
uide losses and inter-resonator coupling strengths, and also 
outperform single ring resonators. We also developed a quan- 
tum mechanical coupled-mode theory which may be appli- 
cable to a generic class of waveguide or resonator based in- 
tegrated photonic quantum light source, and evaluated the 
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expected joint spectral intensities for the apodized and un- 
podized cases. Spectral filtering to isolate individual Bloch 
eigenmodes will help for heralding to a pure state. We also 
introduced a concept of tunability of the output Schimdt num- 
ber, given control over individual inter-resonator coupling co- 
efficients. The special cases of flat transmission and flat group 
delay may give nearly pure heralded states without need for 
additional filtering. 
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